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ABSTRACT 



Context. Starbursts and substantial variations in the star formation histories are a common phenomenon in galaxies. Although predominantly 
found in interacting galaxies, they also occur in isolated galaxies. 

Aims. We study the stability properties of isolated star-forming dwarf galaxies with the aim of identifying starburst modes. The impact of 
the stellar birth function, i.e. a spontaneous and an induced star formation mode, the initial mass function (IMF), the stellar feedback and the 
interstellar medium (ISM) model on the galactic star formation history are investigated. We especially focus on dynamically driven starbursts 
induced by stellar feedback. 

Methods. We apply a one-zone model for a star-gas system coupled by both mass and energy transfer. Additionally, we extend the network 
for active dynamical evolution. This allows for a coupling between the dynamical state of the galaxy and its internal properties, such as star 
formation activity or the thermal state of the ISM. 

Results. While the influence of the dynamics on the total star formation rate is strong, especially with nonlinear stellar birth functions, the 
coupling of the internal properties (gas temperature) on the dynamics is rather limited, because radiative cooling keeps the gas temperature 
well below the virial temperature. Because of short cooling and feedback timescales, the star formation rate is close to the equilibrium star 
formation rates. Quasi-periodic starbursts occur, because star formation follows the variations in the gas density induced by decaying virial 
oscillations. This behaviour is quite insensitive to the nature and the details of the stellar birth description, viz. whether spontaneous or induced 
^ ■ star formation is considered or the IMF is varied. A second type of burst is found as an instability operating when the cooling may drop at very 
' low densities with increasing temperature, in regimes beyond 10 4 K. 

, Conclusions. Bursts of star formation occur during transitory phases, when dynamical equilibrium is established. Then they are quasi-periodic 
on the dynamical timescale. Because of short heating and cooling timescales, the star formation rate follows the equilibrium star formation rate 
corresponding to the actual gas density. 

Key words. Galaxies: evolution - Galaxies: starburst - ISM: general 
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1. Introduction 

Variations in the star formation rate (SFR) in dwarf galaxies 
have been deduced from several observations like direct mea- 
surements of individual star formation histories (Dohm-Palmer 
et al. 2002 ), a large variation in the specific SFR, or short gas 
depletion times (Gallagher & Hunter 1984). Especially inter- 
esting is the group of isolated starburst galaxies, because no 
external triggers can be invoked for their high activity. Van 
Zee (2001) studied a large sample of isolated dwarf irregu- 
lar galaxies. For the majority of these galaxies, she found a 
fairly constant SFR over the past 10 Gyrs indicating that they 
reached a quasi-equilibrium state burning on a long gas deple- 
tion time of about 20 Gyr. However, a small fraction of her sam- 



ple shows high star formation activity that is strongly concen- 
trated in the central region and/or characterised by small radial 
scale lengths, i.e. enhanced gas densities. Already Gallagher & 
Hunter (1984) stressed that the observation of starburst dlrrs 
implies a mechanism to organise large-scale star formation 
which cannot operate all the time due to fast consumption of 
fuel. On the other hand, the instability leading to an episodic 
SFR should only work in a small parameter range, because 
most of the isolated dlrrs seem to evolve in a highly self- 
regulated manner. 

The first theoretical models dealing with star formation 
variations in galaxies were based on closed-box models. For 
reasonable models of the interstellar medium (ISM) it turned 
out that stellar feedback is very effective in suppressing in- 



2 



Ch. Theis and J. Koppen: Starbursts in isolated galaxies. I. The influence of stellar birth function and IMF 



stabilities of the ISM (e.g. Ikeuchi & Tomita [1983] Scalo & 
Struck-Marcell |T986t . These models only allow for a burst- 
like behaviour when a long time-delay (of the order of 10 8 
yrs) between star formation and stellar feedback is introduced. 
Recently Quillen & Bland-Hawthorn (2008 ) emphasised the 
importance of such a long time-delay of the feedback for cre- 
ating episodic star formation in Milky Way type galaxies. The 
drawback of these models is the long time-delay that exceeds 
the short feedback timescales related to massive star forma- 
tion by far. The latter are a few Myr at most (if not quasi- 
instantaneously when taking stellar winds into account). It is 
not clear what has kept the ISM from experiencing feedback 
by stars for such a long time. 

Koppen et al. (1995 hereafter KTH95) investigated a sim- 
ple numerical model motivated by full chemodynamical simu- 
lations. They show that the effective SFR is almost independent 
of the detailed recipe for the stellar birth function, provided a 
negative feedback due to the thermal state of the ISM is con- 
sidered in the stellar birth function. This also holds for a multi- 
phase ISM (Koppen et al. 1998 ). All these box models suffer, 
however, from neglecting galactic dynamics. Additionally, they 
are designed for small galactic regions of a few 100 pc. As a 
result, coherence between different unstable regions and global 
starbursts could not be modelled. 

Recently, several attempts have been made to study the evo- 
lution of isolated dwarf galaxies by taking full stellar and gas 
dynamics, as well as star formation and stellar feedback, into 
account (e.g. Hensler et al. 120041 Pelupessy et al. 120041 I2005T ). 
Though these calculations differ in many details, they agree in 
producing large-scale star formation variations on dynamical 
timescales in some of their models. This indicates the crucial 
role dynamics plays for episodic starburst behaviour. 

Starbursts are often associated with galaxy interactions. 
E.g. most of the ultraluminous infrared galaxies are interpreted 
as the result of galaxy interactions or mergers (e.g. Sanders 
& Mirabel fl 9961 ). It is interesting that the response of galax- 
ies to interaction-induced perturbations is not unique or sim- 
ple. Di Matteo et al. (|2007 ) present a subset of merger events 
created in a numerical survey of galaxy collisions. Their Figs. 
8 and 9 show a variety of reaction patterns of the SFR after 
the first closest approach, ranging from strong starbursts to al- 
most no enhancement in the star formation activity. Another 
example for dwarf galaxy interactions is presented in Bekki 
(2008). Despite this model-to-model variation, statistics of 
large galaxy-pair samples show a clear correlation between in- 
teractions and star formation activity. For example using 12500 
galaxy pairs from the SDSS survey Nikolic et al. (2004) find an 
enhanced star formation rate for projected separations below 50 
kpc and an anti-correlation between the projected separation 
and the star formation rate (for more details, cf. to the review 
by Struck|2006|. 

A major disadvantage of the simulations - either for iso- 
lated or interacting galaxies - is their complexity. This prevents 
a detailed investigation of the parameter space and sometimes 
also hides the key physical mechanisms. Instead of analysing 
dwarf galaxies by such detailed and complex models, we study 
here a set of equations that is as simple as possible that incorpo- 
rates the main generic features of the complex numerical sim- 



ulations. The advantage is a (semi-)analytical treatment and a 
more direct physical interpretation of the implemented physics. 
Moreover, extended parameter studies can be carried out much 
more easily with our set of equations than with the full numer- 
ical modelling. 

This paper introduces an extended version of the analytical 
models of KTH95 that allows us to deal with the dynamical 
evolution of a galaxy. By this, the evolution of the star forma- 
tion activity not only depends on the internal properties of the 
ISM, but also on the dynamical state of the galaxy. As an ap- 
plication of our analytical model, the star formation activity for 
dwarf galaxies in the mass regime of a few 1O 9 M is studied. 

In this paper we investigate the influence on the star forma- 
tion by studying a generalized stellar birth function (composed 
of a spontaneous star formation mode with negative thermal 
feedback and an induced star formation mode), the initial stel- 
lar mass function, and the stellar feedback. 

This paper is organised as follows. In Sect. [2] we introduce 
our set of equations. A reference model and the result of pa- 
rameter variations are presented in Sect. [3] Section |4] contains 
a discussion of the results. 



2. The numerical model 

2.1. The non-dynamical closed box model 

Motivated by the complex and computationally expensive 
chemo-dynamical calculations of Theis et al. (1992) (also used 
in the simulations for dwarf galaxies by Hensler et al. (|2004 )), 
KTH95 analysed a set of closed-box models that was designed 
to be as simple as possible, but still keeps the main qualitative 
features of the chemo-dynamical network (e.g. feedback pro- 
cesses). These simplified equations describe the temporal evo- 
lution of the density of short-lived massive stars (s), of gas (g), 
and of the internal energy (<?) of the gas: 



at t 

at t 
de 

dt 



h(g) S -g 2 MT). 



(1) 

(2) 
(3) 



The gas is consumed by star formation denoted by a stellar 
birth function ^big, T) depending on gas density g and temper- 
ature The temperature is related to the energy density by 
e - bgT with b = 3/2 • k/fj.m p (k is the Boltzmann constant, 
/i the mean molecular weight, and m p the proton mass). The 
second term in (Q]) is a source term for the ISM due to stel- 
lar winds and supernova ejecta replenishing the ISM. Massive 
stars are assumed to return a fraction rj « 90% of their mass 
after a mean stellar lifetime r * 10 Myr. Similarly, Eq. (f2j) de- 
scribes the formation of massive stars, as well as their death 



1 Because the thermal state of all gas components is summarised in 
a single temperature T, this temperature is a "typical" temperature of 
the gas characterising the thermal energy or pressure of the ISM. This 
notion is rather similar e.g. to the notion of the external pressure in 
Blitz & Rosolowsky f2004l ). 
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w 12% is the mass fraction of massive stars derived from the 
initial mass function (IMF)). 

The last equation deals with the energy budget of the ISM 
gas. Its internal energy is characterised by its temperature T . 
Assuming a turbulence-driven ISM it receives heating mainly 
by massive stars (e.g. Ly-continuum radiation or type II su- 
pernovae). The conversion factor h can be constant, if all the 
energy is absorbed, or it may depend on gas density, e.g. for 
radiative heating in an optically thin gas. The main dissipative 
process of the gas is radiative cooling. The energy loss scales 
with g 2 , where the proportionality factor A(T) is the classical 
cooling function. 

The set of equations is closed by a generalized Schmidt-like 
stellar birth function 



x V b (g,T) = V b , sp (g,T) = C n g n -f(T). 



(4) 



This function corresponds to spontaneous star formation. (An 
induced star formation mode will be discussed later in Sect. 
13.31 ) Different to the classical Schmidt-law we allow for a (neg- 
ative) thermal feedback 



f(T) = e 



-T/T, 



(5) 



which reduces the star formation efficiency in case of high gas 
temperatures@ The exponent n is varied between a linear (n = 
1) and a quadratic (n = 2) Schmidt law; the constant C„ is 
chosen to match the star formation rate observed in the solar 
neighborhood. Setting T s to infinity we recover the classical 
Schmidt law depending only on the gas density. 

In case of radiative cooling, thermal equilibrium is quickly 
established due to short heating and cooling timescales 
(KTH95). As a result the effective SFR ^ (which allows for 
a stationary solution of Eqs. (|2]) and (0) is independent of the 
adopted stellar birth function 



¥e(g,T e ) = g 



2 MT e ) 



(6) 



where T e is the equilibrium gas temperature. First, it is remark- 
able that is between a linear and a quadratic Schmidt-law 
depending on h(g), but not depending on especially on the 
exponent n. Secondly, star formation depends almost exclu- 
sively on gas density. This means that such systems undergo 
starbursts only when the gas density is strongly enhanced, e.g. 
by the result of a collapse, an infall, or an induced structural 
change of the galaxy. 

2.2. Including the dynamics 

We extend our closed-box model for dynamical evolution by 
considering a volume with characteristic radius R (e.g. the half- 
mass radius) of the baryonic system of the galaxy. The equation 
of motion - motivated by the corresponding equation for a thin 



2 A negative temperature T s corresponds to a positive feedback. 
However, it is related to the gas temperature (velocity dispersion) 
alone. Therefore, we used a physically better motivated ansatz for pos- 
itive stellar feedback in Sect. 13. 31 



gaseous shell - is given by 
d 2 R 
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(7) 



The first two terms in (O describe the gravity, i.e. the contri- 
bution of the dark matter potential and the self-gravity of the 
baryons with the total mass Mi,. The third term estimates the 
contribution of the pressure P. The fourth term denotes the an- 
gular momentum conservation where j m . dx is the maximum spe- 
cific angular momentum derived for a circular orbit at the initial 
radius R(t — 0) — 7?; n j. The last term describes decaying virial 
oscillations as observed in collapse simulations, and v la d is the 
radial velocity dR/dt. The frictional timescale Tf r i c = Cfri c Tff is 
normalised to the free-fall timescale of the dark matter halo 
at the initial radius. 

From the (half-mass) volume V(t) = 4^/3 ■ R 3 (t) of the 
system and the mass of each component the densities can be 
derived, e.g. the gas density reads 



1 M g (t) 
* (f)= 2'T« 



(8) 



Integrating Eqs. (Q~|) - © over volume yields the equations for 
the total masses of the gas, the massive stars, and the total in- 
ternal energy E = eV inside R: 



dM, 



— = -¥(R,g,T)+^ 

dt t 

dE E dV 

— = h{g) M s P — . 

dt T c00 \(g, T) dt 



(9) 
(10) 
(11) 



The total SFR "V(R,g, T) = y b (g, T) ■ V (in M^r 1 ) de- 
pends not only on the ISM parameters g and T, but now also 
on the current size R of the system. It is worth noting that 
this applies only for a nonlinear star formation law, because 
¥(R,g, T) = VbiMg/V, T) ■ V oc M" g ■ V l -"\ 

The temperature is related to the internal energy E by 
E - 1/2 bM g T . The factor 1/2 stems from the assumption that 
R is the half-mass radius. The last term in Eq. (fTTT > denotes the 
PdV work. The pressure P is given by P = 2/3 E/V — 2/3 bgT. 
In case of radiative dissipation, the corresponding timescale 
T CO oi(^, T) is defined as 



TcoolC?, T) = 



bT 
gMT) 



(12) 



Dissipation by radiative cooling is included by a cooling func- 
tion combined from Dalgarno & McCray ( 1972) for T < 10 4 K 
and Bohringer & Hensler (1989) for higher temperatures (and 
a metallicity Z = l/10Z o ). 

For the dark matter halo, we apply the universal halo sug- 
gested by Burkert ( 1995) for dwarf galaxies. The scale radius 
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ro of the Burkert halo was set to 8 kpc. For a baryonic mass 
of M gas = 2 ■ 1O 9 M and = 8 kpc the initial baryonic mass 
fraction is about 10%. 

The constants Cp, Cj, and Cf,; c are adjusted to a 3d-SPH 
simulation of a collapsing gas sphere in a given static dark mat- 
ter halo potential. The simulations were done with the SPH part 
of the code described in Harfst et al. (2006). Two sets of con- 
stants were derived: set 1 (C P = 1.91, Cj = 0.15, Cf ric = 2.06) 
gives a good fit to the early evolution, whereas set 2 (Cp = 1 .97, 
Cj = 0.11, Cfrk = 0.88) describes the late evolution (Fig. |TJ. 
The main freedom in the choice of parameters comes from the 
damping behaviour that is covered by a variation of a factor of 2 
in the constant Cf,i c , i.e. in the frictional timescale. When min- 
imising the deviations over all periods the best fit was reached 
by C P = 1.74, Cj = 0.13, C fric = 1.64. If not stated differently, 
we apply these values. 

The Eqs. (|7]l - (fT2l yield a closed set of equations includ- 
ing a first approximation for the dynamical evolution. We solve 
the differential equations using a fourth-order Runge-Kutta in- 
tegrator with adaptive timestep (Press et al. 1992). The integra- 
tion accuracy is set to 10~ 10 (or better). For convenience we 
also consider an equation for the mass M r of the long-lived 
low-mass stars and the stellar remnants 



^ = (l-0T(*,*,D + (l-i/)^i. 
at t 



(13) 



The formation of the low-mass stars is described by the first 
term and the production of stellar remnants by the second term 
in Eq. JT3V We neglect the gas return by long-lived low-mass 
stars. This approximation is appropriate for gas-rich systems. 
Therefore, our models cannot provide a realistic evolution of 
the gas mass, if the gas fraction is smaller than a few percent. 
From Eqs. © and (TT0T > we get the solution 



M r (t) = M b - (M g (t) + M s (tj) 



(14) 



describing mass conservation, and Mb is the total initial mass 
in baryons. 

A more accurate treatment would include mass return of 
low-mass stars by planetary nebulae and type la supernovae. 
This would imply the numerically rather expensive computa- 
tion of the instantaneous mass transfer rates from the details 
of the past star formation history for the spectrum of stellar 
masses. For our purpose of studying the general features of 
galaxy evolution with a simplified dynamics, such an elaborate 
approach would not be such an advantage that it would justify 
the computational effort. 

2.3. Units 

For our calculations we used the following units: 1 Mq, lpc, 
and 1 Myr. If not stated differently, parameters and results are 
given in these units or the corresponding derived units, e.g. 
mass densities in Mq pc~ 3 . 
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Fig. 1. Temporal evolution of the half-mass radius in a 3D-SPH 
simulation of a collapsing gas sphere in a dark matter halo 
and the result for the simple equation of motion (0. Shown 
are models fitting the early stage well (with Cp = 1.91, Cj = 
0.15, C f = 2.06, left) or the late stage (with C P = 1.97, Cj = 
0.11,C/ = 0.88, right). 

3. Results 

3.1. The reference model 

As a first example we present a model with a total baryonic 
mass of 2 ■ 10 9 M Q starting at a radius of R m \ - 8kpc. The stellar 
birth function is parametrised by a Schmidt-exponent of n = 
1.5 and a constant Ci.g = 0.06. Additionally, a feedback term 
with T s = 100K is considered. The initial temperature is set to 
10 5 K close to the system's virial temperature 



T = 

J vir — 



Cpk 



G(M DM (R)+±M b ) 



R 



(15) 



This definition is motivated by Eq. ^} assuming force equi- 
librium between thermal pressure at a temperature T v [ r and the 
gravitational forces of a nonrotating system, and Mdm(R) is the 
mass of the dark matter component enclosed within the half- 
mass radius R of the baryonic component. The heating rate is 
set to 10 5I erg per supernova, i.e. h = 2.4 ■ 10 5 pc 2 Myr~ 3 . 

Figure [2] shows the early evolution: Due to strong cool- 
ing, the gas temperature drops almost immediately to values 
near 10 4 K bringing the system far out of virial equilibrium. In 
the ongoing collapse, the density increases. Cooling and den- 
sity enhancement continue, until star formation becomes more 
prominent after about 100 Myr. The first stars reinject energy to 
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Fig. 2. Initial temporal evolution of a model with a baryonic mass of 2 ■ 10 9 M . Shown are the masses of the baryonic components 
(upper left), i.e. gas (solid), low-mass stars, and stellar remnants (dot-dashed), and massive stars (dashed), the mean radius of 
the system (upper right), the total SFR (lower left), and the gas temperatures (lower right). The SFR is compared with the gas 
consumption averaged over 10 Gyrs and the mean star formation up to time ?, i.e. C v P(f')d? / [t. The temperature is compared 
with the virial temperature of the system including the self-gravity of the baryons and the dark matter contribution. 
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Fig. 3. Same as Fig. [2] but for 5 Gyr of evolution. 
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Table 1. Properties of the numerical models 



model 




comment 


A 

A 


i ^ a aa 1 a2 
l.j, U.Uu, 1U 


reference model 


Az 


i ^ o fiA 1 o4 
l.j, U.Uo, 1U 




TWT1 


1 o o 007 ^ 
l.U, U.UU/, co 


— — — - — — — — — — j 

no feedback, Schmidt law 


DiNz 


1.1, U.UU/, co 




DN3 


1.5, U.UU/, co 






9 007 nr. 
Z.U, U.UU / , co 




DIN j 


z.j, U.UU/, co 




LAN D 


9 7 on 

z.u, U. / , co 


— — 

enhanced C 


JJ1 1 


1 n n nm i a5 
l.U, U.UU/, 1U 


with thermal feedback 


DT2 


1 007 1 4 

l.U, U.UU/, 1U 




DT3 


1.0, 0.007, 10 3 




DT4 


1.0, 0.007, 10 2 




DT5 


2.0, 0.55, 10 5 


quadr. Schmidt law with feedback 


DI1 


ind. SF only 


induced star formation mode 


DI2 


model A, T s = 10 2 


combined SF modes 


DI3 


model A, T s = 10 4 


combined SF modes 


DX1 


model A, «f = 6% 


varied constant IMF 


DX2 


model A,£ = 12% 




DWK1 


1.5, 0.06, 10 2 


Weidner-Kroupa IMF, SN heating 


DWK2 


1.5, 0.06, 10 2 


radiative heating 



the ISM leading to the quasi-equilibrium stage known from the 
box models of KTH95. Different to them, the density increases 
because of the continuing collapse: the temperature evolves on 
the line of equilibrium temperatures characteristic of the den- 
sities set by the dynamical state of the system. Though the star 
formation and, by this, the stellar energy injection grow, too, 
the feedback cannot prevent further collapse, because cooling 
is still too efficient to allow the gas to reach the virial tempera- 
ture. The collapse proceeds until t = 0.16 Gyr, when the angu- 
lar momentum conservation in (|7]l becomes dominant and the 
system expands again. By this, the first starburst episode comes 
to an end and the cooling rate drops. The gas is quickly heated 
to 10 4 K; however, the gas never reaches the virial temperature 
thanks to the large increase of the cooling rate at temperatures 
beyond 10 4 K. The next star formation cycle starts when the gas 
is collapsing again. 

The longterm evolution of this model shows that the pre- 
viously discussed bursts are just transitory phenomena caused 
by the initial departure from dynamical equilibrium (Fig. |3j. 
The oscillations are almost completely damped after 2 Gyr. 
However, only the first two or three peaks in the SFR are 
enough strong to be characterised as a starburst by standard 
definitions. Later on no bursts are found. Such behaviour, i.e. 
only an initial burst, was found in almost all models. 

Two general properties of our models can already be seen 
in the reference model. First, despite the coupling of stars and 
gas in the dynamical equation, Eq. ©, the dynamics do (al- 
most) not care about the star formation and the related stellar 
feedback. This is caused by the very efficient cooling resulting 
in gas temperatures well below the virial temperature. Second, 
the SFR is strongly modulated by the dynamics. Exceptions to 
both rules exist, but these are either singular or in physically 
uninteresting regimes (and beyond the validity of our model 
equations). 
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Fig. 4. Temporal evolution of the SFR for models with different 
exponent in the stellar birth function: n — 1 (model DN1, filled 
circles), n = 1.1 (DN2, dot-dashed), n = 1.5 (DN3, solid), 
n — 2 (DN4, open squares) and n — 2.5 (DN5, dashed). The 
normalisation constant C„ is kept constant at 0.007. No nega- 
tive thermal feedback term is included in the stellar birth func- 
tion. 

3.2. Influence of stellar birth function 

In this section we address the influence of the stellar birth func- 
tion and its parametrization according to Eq. @J. In a first se- 
ries of models (DN1 to DN5), we varied the exponent n in the 
stellar birth function fixing the constant C„ to 0.007. We also 
neglected the negative thermal feedback; i.e., we consider the 
limit T s — > oo. This represents the classical Schmidt-law for 
which the stellar birth function only depends on density, but 
not on the temperature of the gas. 

The Schmidt exponent. Figure [4] shows the temporal evolu- 
tion of the SFR for different exponents n. With increasing «, 
the SFR becomes increasingly sensitive to the density vari- 
ations. The peaks of the SFR are fairly similar compared to 
the large SFR variations, whereas the amount of available fuel 
varies strongly with n as a time integration of the SFR shows: 
for the linear Schmidt law, it takes 200 Myr to convert 50% of 
the gas into stars, whereas for n — 1 .5 this takes about 1 .4 Gyr. 
With increasing nonlinearity of the stellar birth function 
the terms containing the gas density become more important 
for the SFR and the coupling to the dynamics is reflected in the 
global SFR. For n — 2 all oscillations of the SFR within the first 
10 9 yr have amplitudes greater than one order of magnitude. 

A linear Schmidt law (n — 1) with no thermal feedback 
term in the stellar birth function represents a singular case with 
respect to the SFR For n = 1, the SFR is proportional to the 
total gas mass, but independent on the gas density: *P = T/, • V — 
C\g ■ V = C\M g . This leads to an exponential decay of the 
gas mass, irrespective of the dynamical state of the system: the 
virial oscillations do not show up in the exponentially decaying 
SFR or the amount of massive stars (Fig. [5] left panels). 

During the early evolution, the gas temperature is always 
close to 10 4 K (Fig. [5] lower right panel). The high tempera- 
tures are a result of stellar heating. Because of a lack of nega- 
tive thermal feedback in model DN1, the temperature does not 




Fig. 5. Temporal evolution of model DN1 characterised by a linear Schmidt law (n = 1, C\ = 0.007), but without a negative 
thermal feedback term in the stellar birth function. For further details see Fig. [2] 
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Fig. 6. Temporal evolution of model DN2 characterised by a Schmidt law (n = 1.1, C„ = 0.007), but without a negative thermal 
feedback term in the stellar birth function. For further details see Fig. [2] 
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Fig. 7. Temporal evolution of model DN4 characterized by a Schmidt law (n = 2, C„ = 0.007), but without a negative thermal 
feedback term in the stellar birth function. For further details see Fig. [2] 



regulate the stellar birth function. Only the strong rise in the 
cooling function beyond 10 4 K prevents the gas from reaching 
higher temperatures. On the other hand, the temperatures are 
too low in the early stage to become dynamically important as 
a comparison with the virial temperatures shows (lower right 
panel in Fig. 0. This lack of star formation induced feedback 
to the dynamics has already been found in the reference model. 
Thus, for the linear Schmidt law, both the dynamics and the star 
formation history are practically decoupled. 
The normalisation constant. In a next step we study the in- 
fluence of the constant C„ in the stellar birth function. For a 
quadratic Schmidt law we varied C„ by two orders of mag- 
nitude. Initially, the difference in the normalisation constants 
is directly reflected in the enhanced SFR of model DN6 (Fig. 
|S}. The enhanced SFR leads to a stronger stellar feedback for 
model DN6. By this, the mean radius increases by about 50% 
and the period of the variations is accordingly longer compared 
to model DN4. The oscillations are more strongly damped and 
already after about 500 Myr the mean SFR is lower for the 
model with the larger star formation constant. The latter is a 
direct consequence of the lack of fuel after the first strong star- 
burst for model DN6. 

The negative thermal feedback. So far, we have neglected 
the thermal feedback exp{-T /T s ) in the stellar birth function 
by setting T s = oo. When we consider high, but finite, temper- 
ature scales T s , the SFR evolves similar to the model without 
any feedback. Figure|9]compares different T s values for a linear 
Schmidt law. In the case of T s = 10 5 K (model DTI) the SFR 
is almost identic to the case of T s = oo (model DN1). After 




0.8 1.2 

time (Gyr) 



2.0 



Fig. 8. Temporal evolution of the SFR for the models DN4 and 
DN6 both characterized by a quadratic Schmidt law with differ- 
ent normalisation constants: C„ = 0.007 (DN4) and C„ = 0.7 
(DN6). 



about 1 .9 Gyr the SFR suddenly ceases caused by the rapid ex- 
pansion of the gaseous system when the density drops below 
a critical density. As discussed in Sect. 14.31 the heating then 
overcomes the high cooling rates beyond 10 4 K and the tem- 
perature rises to values above T s — 10 5 K making the feedback 
term operational. 

Reducing T s to 10 4 K (model DT2) brings T s closer to 
the actual temperatures of the ISM mainly set by the strong 
rise in the cooling function. Variations in the temperature in- 
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Fig. 9. Temporal evolution of the SFR for models with a lin- 
ear Schmidt law (n = 1) and a negative thermal feedback term 
with different values of the temperature parameter T s : T s = oo 
(model DN1, no thermal feedback in the stellar birth function, 
filled circles), T s = 10 2 K (DT4, solid), T s = 10 3 K (DT3, 
dot-dashed), T s = 10 4 K (DT2, short-dashed), and T s = 10 5 K 
(DTI, long-dashed). The normalisation constant C„ is kept con- 
stant at 0.007. 



fluence the stellar birth function much more strongly than in 
model DTI . This can be seen in the reduced overall SFR, which 
shows a modulated exponential decay. Setting T s to 10 2 or 10 3 
K makes the density variations and its impact on the tempera- 
ture by cooling to the dominant factor for the global SFR. The 
oscillation periods are then given by the dynamical timescale 
and the star formation reaches a quasi-equilibrium state when 
the virial oscillations have been decayed. 

It is also interesting to consider the long-term evolution of 
model DT2 (T s = 10 4 K, Fig. OH: after about 6 Gyr of self- 
regulated quasi-equilibrium evolution, the SFR enters an oscil- 
lating period over the next 3 Gyr. Then again the evolution be- 
comes rather stable. This variation in the SFR is directly caused 
by the temperature fluctuations between about 1.4 ■ 10 4 K and 
8 • 10 4 K as the upper two panels in Fig. [TOl show. Since T s 
amounts to 10 4 K, the feedback term varies by about a factor 
of 10 3 . A comparison of the radial variations and the tempera- 
ture variations shows that the fiuctations are not caused by the 
dynamical evolution (Fig. [10] lower panel): the radial oscilla- 
tions are too small in amplitude and are not in phase with the 
temperature variations. 

Though e.g. the temperature variations are too large to al- 
low for a complete explanation by linear perturbation theory, 
a linear stability analysis might give a first clue to the exis- 
tence of this unstable regime. Since the dynamical evolution is 
unimportant here, we can apply the results of the stability anal- 
ysis for the nondynamical set of equations ([]]) - (0) given in 
KTH95: a necessary condition for stability is that Eq. (28) of 
KTH95 holds. In the limit of low gas densities, we can rewrite 
this stability criterion as 

T C ooi>-A'(r)r. (16) 



This is the case, if the logarithmic slope A' of the cooling 
function is positive definite (as assumed for the whole tem- 
perature regime considered by KTH95). However, the more 
realistic cooling function applied here drops for temperatures 
beyond the peaks caused by hydrogen and helium line emis- 
sion. Though this behaviour depends on metallicity (becoming 
weaker with increasing metal content), it persists for metallic - 
ities ranging from metal-free to solar metallicity (Bohringer & 
Hensler 1989). Thus, the stability of the system is not guaran- 
teed anymore. 

Assuming that the logarithmic slope of the cooling func- 
tion is close to one, instability can only occur if the cooling 
timescale is longer than the mean stellar lifetime r of the mas- 
sive stars. In the early stages when the gas density is high, the 
instability criterion will not be met due to the high efficiency 
of cooling. However, in the late stages when most gas is con- 
sumed (or lost), the cooling timescale becomes long even for 
temperatures above 10 4 K. 

For a metallicity of 1/10 Zq, the cooling curve according 
to Bohringer & Hensler has a first region of negative slope be- 
tween 1.5 TO 4 K and4T0 4 K. Therefore, the evolution becomes 
unstable once the equilibrium temperature rises to 1.5 • 10 4 K. 
The next unstable region caused by He line emission exists for 
equilibrium temperatures close to 10 5 K. However, in model 
DT2 this regime is not reached within the first 10 Gyr. 

3.3. Influence of positive stellar feedback 

The star formation recipe discussed in Sect. 13.21 corresponds 
to a spontaneous star formation mode with a negative stellar 
feedback. However, there might also be positive stellar feed- 
back (see e.g. the stochastically self-propagating star formation 
models (SSPSF) by Gerola & Seiden [1978j >. A model often dis- 
cussed is the SN-induced star formation (e.g. Ehlerova et al. 
[19971 Hosokawa & Inutsuka l2"005l 12006b : the basic idea is that 
expanding SN-shells sweep up the ambient ISM in a thin, very 
dense shell. Such shells cool very efficiently which results in 
fragmentation and subsequent star formation. 

To investigate the influence of an induced star formation 
mode, we extend our stellar birth function to 

V b (g, T; s,R) = ^, sp fc T) + W bM (g, s,R) , (17) 

which allows for both a spontaneous star formation mode and 
an induced star formation mode. For the spontaneous star for- 
mation mode ^b^pig, T) we adopt the description applied in the 
previous sections; i.e., Eq. ©. For the induced star formation 
rate we estimate the gas mass swept up by SN shells and turned 
into stars by 

TynCfc = — ■ MR^s,g),R) , (18) 

where /■ denotes the fraction of a galaxy covered by SN shells, 
i.e. gfi is the maximum amount of gas undergoing SN shell- 
induced star formation. A fraction rj, of this mass is then con- 
verted on a timescale t,- into stars. The factor is related to 
the volume R 3 , filled by massive superbubbles undergoing frag- 
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Fig. 10. Late temporal evolution of the gas temperature (upper 
and lower panel, solid line), the SFR (middle panel) and the 
mean radius of the system (lower panel, dashed line) for the 
model with a linear Schmidt law (n = 1) and a negative thermal 
feedback term with T s (DT2). The lower panel shows only a 
short fraction of the time displayed in the upper panels. 



mentation relative to the volume of the galaxy. For simplicity 
we assume 



MR*(s,g),K)=l-e 



-0W«r 



(19) 



In the case of low star formation activity, only a small fraction 
of the galactic volume is affected by induced star formation, 
whereas for high star formation activity most of the gas mass 
might be swept up. The related star formation timescale t,- is 
given by the time for bubble creation. The efficieny factor 77, 
gives the fraction of gas accumulated in shells that is converted 
into stars. We used a value of 77, = 0.1 here. The size and the 
timescale for the bubble evolution were taken from Eqs. (13) 
and (14) of Ehlerova et al. (119971) . It should be noted that this 
star formation recipe is only meant to be a simple qualitative 
estimate, not a detailed sophisticated induced star formation 
model. Anyway, we believe that this ansatz covers all essential 
ingredients for the SN-induced star formation mode. 
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Fig. 11. Temporal evolution of the SFR for models with dif- 
ferent stellar birth modes: induced star formation only (dashed 
line, DI1), spontaneous star formation only (solid line; T s - 
100 K (model A, upper panel) and T s = 10 4 K (A2, lower 
panel)) and both star formation modes together (T s = 100 K 
(DI2, plus signs in upper panel) and T s = 10 4 K (DI3, open 
triangles in lower panel)). 
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Fig. 12. Temporal evolution of the relative contribution of the 
spontaneous (solid) and the induced (dot-dashed) star forma- 
tion mode to the total mass converted into stars for model DI3. 
For comparison the combined SFR is also shown (plus). 
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Figure QTJ displays the SFR for different combinations of 
the star formation modes. In case of a low temperature T s = 
100 K, the overall SFR is dominated by the induced star for- 
mation (ISF) mode (upper panel of Fig. ITTb as the agreement 
between model DI1 (ISF only mode) and the model DI2 (com- 
bined star formation modes) demonstrates. This is caused by 
the high temperatures of about 10 4 K created by the massive 
stars formed by induced star formation. Due to negative ther- 
mal feedback the spontaneous star formation is almost com- 
pletely suppressed. This suppression can, of course, not take 
place in a purely spontaneous star formation model (model A). 
In that case, the amplitudes of the density driven SFR varia- 
tions increase because of the temperature dependent variations 
in the feedback term. The dominance of the dynamical evo- 
lution is reflected in the almost identical periods of the SFR 
in all models. The slightly longer periods for the models DI1 
and DI2 (which include induced star formation) stem from the 
larger equilibrium radii leading to a longer free-fall time. 

To have a substantial contribution from the spontaneous 
star formation mode in a stellar birth function including in- 
duced star formation, the high efficiency of the thermal feed- 
back has to be reduced. This can be realized by an increase in 
the feedback temperature scale T s close to or above the actual 
gas temperature. Model DI3 has a T s of 10 4 K: models using ei- 
ther of the two star formation modes (induced SF, spontaneous 
SF) or both of them become then comparable in star formation 
amplitudes, peak values and the overall production rate of stars 
(Fig. QT| lower panel). Still, the SFR is dominated by the dy- 
namically driven oscillations. It is interesting to note that after 
1 Gyr the combined SFR (model DI3) becomes smaller than 
both single-mode models (DI1 and A2). 

The relative fraction of the spontaneous and the induced 
star formation mode in the combined model DI3 with respect 
to the total number of stars formed is shown in Fig.Q~2] at the 
very beginning the spontaneous star formation exceeds the in- 
duced star formation. However, with an increasing overall SFR 
the induced star formation begins to dominate and the spon- 
taneous star formation drops to 30%. The contribution of the 
induced star formation drops already before reaching the first 
peak of the SFR: this is caused by the dynamics, i.e. the in- 
crease in the gas density during the collapse of the system and 
the subsequent enhanced spontaneous star formation. The tem- 
perature of the gas remains almost constant close to 10 4 K dur- 
ing this evolution. After a Hubble time about 62% of the stars 
are produced by induced star formation. 

It is remarkable that adding induced star formation does 
not (necessarily) result in a quick consumption of all avail- 
able gas. For our fiducial model, roughly half of the stars are 
formed in one of the two star formation modes. Though this 
fraction can be varied by e.g. amplifying the thermal feedback 
term of the spontaneous star formation mode or changing the 
efficiency factor 77, of the induced star formation, none of the 
star formation modes leads to a completely different behaviour 
of the system: still the overall star formation rate is controlled 
by the global dynamics. 
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Fig. 13. Temporal evolution of the SFR for the reference model 
A (solid), models DX1 and DX2 with a constant IMF-derived 
fraction £ of massive stars (£ = 0.06: open squares; £ = 0.24: 
filled circles) and the Weidner-Kroupa type variable IMF mod- 
els based on SNII heating (DWK1, dashed) and radiative heat- 
ing (DWK2, dot-dashed). 

3.4. Influence of the IMF 

Constant IMF. A variation in the stellar initial mass function 
will mainly affect the mass fraction £ of massive stars. To test 
limiting cases we varied the time-independent £ by a factor of 
2. Though such a change is practically beyond the limits of 
generally adopted standard IMFs, the impact on the evolution 
is small; e.g. the SFR is qualitatively identical to the reference 
model for both ^-values (Fig.[T3l. The small quantitative differ- 
ences (e.g. lower SFR for larger £) can be understood as a result 
of enhanced (^ = 0.24) or reduced thermal feedback. The weak 
response of the SFR (and especially its peak values) is a result 
of the strong self-regulated coupling of the star-gas system. As 
a result, the evolution of the star-gas system is practically in- 
sensitive to reasonable variations in the constant IMF. 
Variable IMF. Recently, Weidner & Kroupa (120051 ) suggested 
a variable IMF that depends on the global star formation activ- 
ity of the galaxy. Their main idea is that relations exist between 
the global galactic SFR and the maximum mass of the molecu- 
lar clouds as well as between the cloud mass and the maximum 
mass of the most massive stars in star formation regions. By 
this, the global SFR affects the IMF for massive stars which 
are responsible for the main heating and the metal production. 

To investigate such a scenario we modified our equations 
by introducing an efficiency factor /wkOP) for the production 
of massive stars; i.e., we replaced £ by ^•/wk( 1 P)- The efficiency 
factor was given by a simple fitting formula adopted from the 
analysis of Koppen et al. (120061) : 

I 1 - 0.8 • e~ xl2 x > 

/wk = | for (20) 

{ 0.2 e x x<0 

with x = 3 + log [*F/(M yr -1 )]. One should note that the effi- 
ciency factor only varies slightly over the physically interesting 
regime of x > or *F > 10~ 3 M yr~' . 

We considered two cases for the energy feedback. In the 
first case we assumed that the feedback is given by type II su- 
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pernovae (model DWK1). In that case the heating coefficient 
corresponding to the energy input per supernova has to be mul- 
tiplied by the number of massive stars. The number of massive 
stars, however, depends mainly on the lower mass limit of the 
mass range of massive stars, and not on the upper mass region 
that is mainly affected by a Weidner-Kroupa type IMF. A de- 
tailed calculation shows that the conversion factor between the 
number of massive stars and their mass varies only by a factor 
of 2 when changing the SFR from 1CT 3 M yr~' to 100 M Q yr~' . 
In our model DWK1 we neglect this small variation and keep 
the energy equation. 

In the second case, we assume that the heating is mainly 
done via radiation (model DWK2). Then, the heating coeffi- 
cient h becomes a function of the SFR, because of the strong 
dependence of the stellar radiative energy input on the upper 
stellar masses caused by the mass-luminosity relation. A de- 
tailed calculation shows that the corresponding efficiency fac- 
tor is similar to the one for the mass fraction. Therefore, we 
use for simplicity the same factor /wk, i-e. h is replaced by 
h ■ /wkC^) in Eq. (flTT i (or Eq. J3), respectively). 

Figure Q~3] shows the evolution of both models. The dif- 
ference between models with and without a variable IMF are 
marginal in both, amplitude and timing. Compared to the ref- 
erence model A, the Weidner-Kroupa-type IMF models have 
a slightly enhanced SFR. This enhancement is necessary to 
balance the deficit in the stellar feedback caused by the re- 
duced number of massive stars in the Weidner-Kroupa-type 
IMF models. 

In general the qualitative evolution of the models is not af- 
fected by the shown IMF variations, either if we change the 
mass fraction £ assuming a stationary IMF or if we apply the 
temporal variable star formation-dependent IMF suggested by 
Weidner & Kroupa (2005). The main reason is the very efficient 
self-regulation. 

4. Discussion 

4.1. Timescales 

Basically four physically distinct timescales are involved in our 
description: the heating and cooling timescales, ta, and Tdiss, 
the dynamical timescale Td yn and the gas consumption (or star 
formation) timescale tsf. KTH95 already has shown that tsf 
is usually longer than the timescales governing the energetics 
of the ISM. This also holds in case of a multiphase ISM model 
(Koppen et al. I1998I I, where the timescale for the additional 
equilibrium between evaporation and condensation of the gas 
clouds is well separated from the other two, leaving them un- 
perturbed. Therefore, we expect that inclusion of a multiphase 
ISM model would not alter the basic behaviour of our system. 

For typical mean gas densities the dynamical timescale ex- 
ceeds also the cooling and heating time (though it is less than 
the gas consumption time). Therefore, the system almost in- 
stantaneously adjusts to the equilibrium state given by the ac- 
tual gas density. Provided the temperature remains in a regime 
with a positive slope of the cooling function, the system is then 
stable. This holds e.g. for temperatures below 10 4 K, (which is 
practically the temperature regime for high and moderate gas 
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Fig. 14. Evolution of the reference model A (solid) and the 
model DT5 (quadratic Schmidt law with feedback term, T s = 
10 5 K; dashed line) in the phase space R-v ra i. The inset shows a 
blow-up of model A near its equilibrium position reached after 
the initial oscillations. 
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Fig. 15. Evolution of the reference model A in the R-T plane. 
The inset shows a blow-up near the equilibrium position 
reached after the initial oscillations. In the inset linear temper- 
ature values are used. 



densities). Thus, the radial oscillations are only initial transient 
virial oscillations acting on a dynamical timescale. 

This behaviour turned out to be rather robust against vari- 
ations in model parameters like stellar heating, gas mass, ini- 
tial spatial extension, stellar birth function (with the exception 
of a linear stellar birth function resulting in an exponentially 
decaying SFR). Also different IMF parametrizations or even 
a variable IMF (like the integrated galactic IMF suggested by 
Weidner & Kroupa 2005) do not change the behaviour of the 
system qualitatively provided the feedback mechanisms are not 
cutoff. 

4.2. Generic behaviour 

Our calculations show that the evolution of the coupled star- 
gas system is dominated by the dynamical state of the system. 
Due to the short heating and cooling timescales, the ISM re- 
laxes quickly. Thus, the system's behaviour is governed by the 
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dynamical equation Eq. (0, which has the structural form 
dR 



Tt = Vrad 



dv 



rad 



dt 



A 2 Ay A 3 

7 + + — 7 - A 4 v i 

R 2 R R 3 



(21) 
(22) 



with coefficients A, that are positive for physically sensible de- 
scriptions. For simplicity we only consider the self-gravity term 
here. With a proper choice of the time unit, A4 can be set to 
unity. 

For v ra d = there exist equilibrium positions, where the 
forces cancel out 



Ai A 2 A 3 
— - ^7 + —t = 



Re Ri 



or 



Ri 



A x R 2 e - A 2 R e + Ai = 
with the formal solution 

Re± = ( A 2 ± > /a2-4AiA 3 ) . 



(23) 



(24) 



(25) 



If A\ > 4A1A3, two real equilibrium radii exist. Note that for 
Ai = (i.e. a gas temperature well below the virial tempera- 
ture) one simply has a single solution 

A3 



R, = 



(26) 



The analysis of the local stability of these equilibrium 
points (via the linearized equations) of the Jacobi matrix J 



dR 
dR 



= 1 



<9v ra( j 

dVrad = _M 2A 2 3A 3 

dR Rl + Rl R A e 

3v la d 



dv md 
with 



-1 



Trace J 
detj 

Discr = 1 -4detJ/(TraceJ) 2 



-1 

A, 2A 2 
Rl~~kT 



3A 3 
Ri 



(27) 
(28) 
(29) 
(30) 

(31) 
(32) 

(33) 



shows if equilibrium solutions exist (A\ > 4A1A3) that R e - al- 
ways is an attracting focus, since the trace and the discriminant 
are negative for all positive coefficients At. Likewise, the other 
point R e+ always is a saddle point (negative determinant). 

The consequences of this general structure are shown by a 
few representative trajectories in the phase space. In Fig. [14] we 
depict these already for the complete system of our equations. 
One notes that if the initial conditions are close enough to the 
attracting focus, the system performs damped nonlinear oscil- 
lations before settling in the equilibrium. If the initial radius or 



the initial speed is too large, there will be a collapse, followed 
by an expansion. This takes the system to the vicinity of the 
saddle point, whence the expansion continues without bounds. 

If one looks more closely at the behaviour in the equilib- 
rium state (shown in the inset in Fig. fT4l) . one finds that the 
equilibrium radius does not stay constant, but always increases, 
more or less rapidly depending on the coefficients A,. The rea- 
son is that the coefficient A 1 of the pressure term in the full set 
of equations depends on the gas temperature. When the sys- 
tem reaches the (dynamical) equilibrium, the gas continues to 
be consumed by star formation. The cooling of the gas there- 
fore slows down, and the gas temperature increases slowly, giv- 
ing rise to a slow expansion that contributes to the decrease in 
the gas density. This rise in temperature causes an increase in 
the equilibrium radius, which the system simply follows. This 
leads to the slow drift seen in the inset of Fig. FBI 

In Fig. Q3] we plot a trajectory in the radius-temperature 
plane. One notes that the oscillations in radius and temperature 
go with a nearly constant phase lag, but once the oscillations 
have damped out, a steady secular increase in both radius and 
temperature takes over (cf. to inset in Fi g . [Tsb . 

4.3. Heating-cooling equilibrium at low gas densities 

The equilibrium between heating and cooling plays a crucial 
role for the overall evolution of our system. Cooling is very effi- 
cient for typical galactic mean densities, thereby compensating 
for the stellar heating processes. However, when the gas den- 
sity drops, the heating by stars (hM s oc WjV oc M g ) can over- 
come the cooling (Ag 2 V oc gM g ), so the gas quickly heats up 
(here we assume a linear stellar birth function). Once temper- 
atures become higher than the virial temperature, the gaseous 
system expands and a mass loss occurs even without any dra- 
matic event like a starburst. The critical density g] oss for gas loss 
can be estimated from the balance between heating and cooling 
from Eq. © with s « ?¥t = (C n g" e -777 "' 



5 loSS 



h^C„re 



-TIT, 



A(D 



(34) 



For the values of model DN1 (linear Schmidt law), we 
get gloss ~ 8 • 10~ 5 Mqpc~ 3 by adopting a cooling rate 
of A KT 2L5 erg cm 3 s~' characteristic of temperatures well 
above the drop near 10 4 K. The g\ oss becomes smaller with in- 
creasing Schmidt exponent n. Practically, such low gas densi- 
ties might not be relevant, e.g., because of the gas replenish- 
ment by low-mass stars. 

It is interesting to note that no critical density exists for a 
quadratic Schmidt law. If n is higher than 2, the system would 
even run into a cooling catastrophe. In the case of optically thin 
heating, h in Eq. ( f34b has to be replaced by h g. Thus, the criti- 
cal exponent for a cooling catastrophe already becomes n = 1. 
This means that, for all reasonable stellar birth functions, the 
gas cools down below a critical density. In that case other heat- 
ing mechanisms providing a steady energy source independent 
of the amount of massive stars (hence, the SFR and the gas 
density) might become significant. 
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Fig. 16. Temporal evolution of the mean radius of model DN1 
(solid) and the same model, but neglecting the pressure term in 
the dynamical equation Eq. (dashed). 



4.4. Are the dynamics influenced by stellar feedback? 

In this paragraph we discuss the impact of stellar feedback 
on the dynamical evolution. The latter is described by Eq. (0 
which contains a pressure gradient term: l/g ■ dP/dr ~ T/R. 
Stellar feedback can influence the dynamical evolution via this 
term by heating. Practically speaking, this happens when the 
pressure term is close to or exceeds the gravitational accelera- 
tion terms (the first two terms in Eq. 01), i.e. when the gas tem- 
perature is close to the virial temperature or higher. However, 
in case of the evolutionary stages shown in this paper the gas 
temperatures are always of the order 10 4 K or below, whereas 
the virial temperature is at least a few 10 5 K. Physically, the 
low gas temperatures are caused by the steep rise in the cooling 
function beyond 10 4 K (cf. e.g. Fig. 1 in Bohringer & Hensler 
1989). For typical gas densities, unrealistically high heating 
rates are required to reach (and keep) temperatures beyond 10 4 
K. 

Therefore, the dynamics are almost unaffected by the ther- 
mal state of the gas. This is nicely illustrated in Fig. [16] where 
the mean radius of model DN1 is compared with a model ne- 
glecting the pressure term of Eq. 01: no substantial difference 
between both models is discernible. Small differences develop 
at maximum compression when the gravitational terms and the 
angular momentum almost cancel. For this short moment, the 
small pressure-gradient term contributes recognisably. It acts as 
a reduction of the gravity, e.g. like a density reduction, which 
results in slightly longer periods. 

The situation only changes when either the gas tempera- 
ture rises or the virial temperature decreases. The latter can 
happen for systems of lower mass, e.g., for compact systems 
like globular clusters (no DM, ~ 10 6 Mq) or more extended 
systems like dSph galaxies, when the virial temperatures drop 
down to a few 10 4 K (we will discuss this in a forthcoming 
paper). Alternatively, the gas temperature might be heated to a 
few 10 5 K. Practically speaking, this is only the case where the 
criterion for run-away heating, Eq. Q4I) . is fulfilled, i.e. when 
the gas density drops below a critical value due to gas con- 



sumption. The actual critical values, however, are too low, to 
be practically relevant within the limits of our model. 

The physical reason for the general insensitivity of the dy- 
namics to the thermal feedback induced by stars is the steep rise 
in the cooling function near 10 4 K. As long as radiative cooling 
is the main agent of energy dissipation and as long as the gas 
temperature is kept below 10 4 K, thermal feedback will only 
affect the dynamics of systems with a virial temperature near 
10 4 K. However, for short periods (especially at early stages), 
the gas temperature might reach or exceed temperatures of 10 5 
K briefly affecting the dynamics of systems with virial temper- 
atures of that order. 

However, thermal pressure (and therefore classical cooling) 
does not need to be the main agent counter acting gravity. If 
the system is composed of sufficiently compact objects, their 
velocity dispersion will provide a kinetic pressure, e.g., in case 
of a cloud system, the internal temperature of the clouds is ir- 
relevant for their overall dynamics, which is driven by the ve- 
locity dispersion. In that case, inelastic cloud-cloud collisions 
provide the relevant energy dissipation mechanism. Different 
from radiative cooling, there is no steep increase in the cloud- 
cloud "cooling" function and the arguments preventing a strong 
mutual coupling between dynamics, star formation, and stellar 
feedback are not applicable. Additionally, the timescale for en- 
ergy dissipation due to cloud-cloud collisions is much longer 
than the radiative cooling timescales, which also might allow 
for more of a variety of responses by the system. We will dis- 
cuss the related aspects in a future paper. 

4.5. Comparison with "full" dwarf galaxy simulations 

How do our simple models compare with detailed numerical 
simulations? Though important, it is a difficult question. At the 
moment there are only a small number of sophisticated simu- 
lations available. They differ by many aspects between them- 
selves and from our calculations: e.g. with respect to the imple- 
mented physics, the numerical methods or the initial conditions 
and the general physical set-up. Unfortunately, a standard test 
for dwarf galaxy evolution/physics does not exist. 

With all these differences in mind a detailed comparison 
is not meaningful, but a qualitative comparison should shed 
light on the applicability of our ansatz; e.g., Pelupessy et al. 
( 2004, hereafter PWI04) model the evolution of a rotating disk- 
like dwarf galaxy embedded in a dark matter halo. They apply 
a sophisticated ISM model including heating, cooling and an 
ionisation/recombination scheme. Their SFR develops a quasi- 
periodic behaviour with a period of about twice the dynamical 
time for (weak) bursts. PWI04 interpret this behaviour as the 
result of a kick that drives the gas out of the disk and the sub- 
sequent fallback of the gas initiating the next star formation 
episode. This concept is very similar to our adopted coupling 
between feedback and dynamics, though we are investigating a 
spherical system. 

The overall star formation peaks in PWI04 are much 
smaller than in our calculations (by one order for a fiducial 
model including a feedback term). This might be caused by 
the different initial set-ups, because the spherical configura- 
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tion we study allows for a much deeper collapse than a rota- 
tionally supported disk. As a result, the density and, thus, the 
star formation variations are much higher in our case than for 
PWI04. For comparison, the full Id chemo-dynamical mod- 
els of Hensler et al. (120041 afterwards HTG04) model spher- 
ical galaxies (but without dark matter) and they exhibit SFR 
comparable to ours. Similar to our calculations, the models by 
HTG04 show a strong trend toward self -regulation that has also 
been found in the 3D chemo-dynamical N-body-SPH models 
of Berczik et al. d20021 l including dark matter. The PWI04 mod- 
els exhibit no strong starbursts as the moderate variation in their 
star formation rates shows (Figs. 4 and 6 in PWI04). In contrast 
to them we were looking for strong variations that are rarer, but 
that exist (van Zee 1200 II) . 

Another interesting example of SFR variability is the SPH 
simulation of Struck (120051 hereafter S05), who investigated 
the evolution of a disk galaxy of about 2 ■ 10 1 1 M Q . The gaseous 
component in this simulation was subject to star formation, 
stellar feedback (by SNII) and dissipation (radiative cooling). 
S05 finds a strongly variable star formation with two modes: 
local, off-centre gas concentrations can lead to local, "prema- 
ture" bursts preventing the system from a global starburst. The 
result is an incoherent star formation pattern with an irregular, 
mildly variable SFR. Sometimes, however, large-scale motions 
coherently funnel matter into the central region where a strong 
starburst is then ignited (cf. Figs. 1 and Al in S05). The re- 
lated timescale is about twice the free-fall timescale of a par- 
ticle at 2-3 kpc above the plane. This burst mode corresponds 
physically to the oscillatory "breathing" mode seen in our cal- 
culations, though the masses of the models differ by one order 
of magnitude. Similarly, the burst duration (about 10 Myr) is 
comparable to the burst duration in our models. On the other 
hand, our simple model is unable to resolve purely local effects 
leading to incoherent star formation modes. 

Though a detailed comparison between our simple model 
and the different complex numerical simulations is unreason- 
able because of differences in implemented physics and in ini- 
tial conditions, the qualitative agreement is obvious: the bursts 
essentially stem from the oscillations of the global dynamics. 
The gas simply has to follow them, so the starbursts could ap- 
pear in any system. The tight coupling between the large-scale 
dynamics and the occurrence of global (strong) starbursts is 
corroborated in all the different studies, by this motivating and 
strengthening our ansatz. 

4.6. A speculation about bursts in interacting 
galaxies... 

By definition our model refers to isolated systems. However, 
the initial conditions leading to a collapse of the system could 
also stem from galaxy interactions. Therefore, one might spec- 
ulate about the possible implications for interacting galaxies. 

An important result of our simulations is the strong cou- 
pling of the SFR to the dynamics (with the discussed exception 
of linear Schmidt law without feedback term). To create a star- 
burst, the gas density must be strongly enhanced, otherwise the 
self-regulation prevents large SFR variations. The availability 



of large amounts of gas alone is not sufficient for a starburst as 
the example of poststarburst E+A galaxies with a large amount 
of HI gas demonstrated (Buyle et al. 2006). 

For interacting galaxies, one might speculate that only 
those perturbations leading to structural rearrangements, e.g. 
a substantial inflow of gas or a merger, end up in a star- 
burst. Whether galaxies produce sufficient density enhance- 
ments, hence a starburst, depends not only on the existence of 
a proper perturbation (like an interaction or an internal bar) but 
also on the internal dynamical stability of the involved galaxies. 
This might explain the wide variety of star formation responses 
found in simulations of interacting galaxies (e.g. in Di Matteo 
etal. l2007l ). 

4.7. Simplifications of the model 

Describing galaxy evolution by a few equations always requires 
simplifications. Here we want to comment on a few simplifica- 
tions (or neglections). 

Dynamics. The dynamics of the system is summarily described 
by a mean radius and its temporal change. A decomposi- 
tion into several structural components (like a disk-bulge-halo 
structure) is beyond the limits of this model. Moreover, density 
gradients or even an asphericity are not considered. An exten- 
sion of our simple model might be possible. However, with re- 
spect to structure formation, we think that the related physics 
are too complex to be condensed reliably into another "simple" 
equation. 

Implicitly we also assumed that the dynamics of the differ- 
ent baryonic components can be modelled by the same mean 
radius. This need not to be true, of course. Strictly speaking 
this is only valid when the dynamics of stars and gas is ei- 
ther subject to the same acceleration terms (e.g. similar pres- 
sure/velocity dispersion terms) and/or when the self-gravity 
of the stars is not important for the dynamical evolution of 
the gas. The latter is the case in gas-rich galaxies or in dark- 
matter-dominated potentials (when the self-gravity term re- 
mains unimportant). Since we are interested here in starburst 
galaxies that are mainly gas-rich dwarf galaxies, the restric- 
tion to gas-dominated stages seems reasonable. The dynami- 
cal equations can be easily extended to several distinct compo- 
nents, if required. 

ISM description. In our approach we adopted a single ISM 
phase. An extension to a multi-phase ISM model like in 
Koppen et al. d 19981 ) would be possible and is envisaged. 
However, we do not expect significant changes, because the 
timescales related to the additional multi-phase ISM processes 
(like condensation, evaporation) are separated well from the 
other timescales: they are much longer than the dominant heat- 
ing and cooling processes, so that thermal equilibrium is es- 
tablished rapidly and independently of the fractional ratio of 
clouds and intercloud gas. On the other hand, evaporation and 
condensation of clouds occur much faster than the gas con- 
sumption due to star formation. Thus, we expect that the sys- 
tem would again break down in a hierarchy of well-separated 
equilibria, and the related equilibrium SFR would thus still be 
governed by the equilibrium temperature in the cloud phase. 
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A new effect that arises for a clumpy ISM is another en- 
ergy dissipation process. In the case of a diffuse ISM, dissipa- 
tion is mainly driven by radiative cooling; however, in case of 
a system composed of molecular clouds, inelastic cloud col- 
lisions provide a "cooling" mechanism dissipating the kinetic 
energy related to the motion of the clouds. Though the scal- 
ing of this dissipation process is identical to radiative cooling 
(i.e. oc g 2 ), the related timescales differ strongly. In general, 
cloud-cloud collisions operate on a much longer timescale than 
radiative cooling or even the dynamics. Thus, the dynamical 
evolution of a cloud system might operate not only on a dy- 
namical timescale, but also on a dissipational timescale which 
then directly affects the SFR. We will discuss the implications 
in a forthcoming paper in detail. 

Heating. For the heating we only considered massive stars and 
instantaneous feedback. Though massive stars are considered 
to be a major source of ISM heating, e.g. with respect to tur- 
bulence, other sources might be interesting, too. For example, 
type I supernovae provide a nonlocal and a non-instantaneous 
energy source. Especially, the delayed energy feedback might 
be an important energy source after star formation has ceased. 
Similarly, the energy injected by an AGN might be an impor- 
tant energy source, too. 

5. Summary and conclusions 

We investigated the evolution of star-forming dwarf galaxies by 
means of an extended one-zone model. In contrast to previous 
one-zone models, the dynamics of the galaxy has been taken 
into account by an additional equation of motion for the mean 
radius of the galaxy. Comparison with SPH models shows that 
this rather approximate treatment is reasonable and provides a 
sensible description. This approach allows for a coupling be- 
tween the dynamical state of the galaxy and its internal proper- 
ties, such as star formation activity and the thermal state of the 
interstellar gas. Here we focussed especially on the conditions 
under which starburst episodes could occur in isolated galaxies. 

We find that the seemingly complicated system of equa- 
tions for the numerous physical variables (radius, radial speed, 
gas mass, star mass, thermal energy) breaks down into quite 
a clearly separated system for the internal conditions of the 
galaxy and the dynamical aspects that evolve more slowly. 
Thus, star formation and other internal processes are strongly 
governed by the dynamical state. In particular, the total SFR 
follows the evolution of the volume directly, if it depends non- 
linearly on the local gas density. On the other hand, the dy- 
namics are only very weakly influenced by the internal state, 
because the gas temperature remains well below virial temper- 
ature, primarily because of the strong rise of the radiative cool- 
ing function near 10 4 K. The consequence is that the evolu- 
tion of the dynamical state of the galaxy remains fairly robust 
against changes and details in the internal condition, specifi- 
cally the details of the prescriptions for star formation. We also 
expect that the inclusion of a multiphase description of the ISM 
would not upset this characteristic. 

The evolution of the total SFR is thus dominated by the 
damped virial oscillations that cause enhancements of the SFR 
whenever the galaxy is in a more compressed state. This main 



type of quasi-periodic starbursts is related to the dynamical 
timescale of the system. Typically, the star formation varia- 
tions follow the variations in the gas density induced by decay- 
ing virial oscillations. Because of the short heating and cooling 
timescales the established SFR remains close to the equilibrium 
SFR determined by the current mean gas density. Modifications 
in the adopted star formation modes, i.e. assuming both sponta- 
neous and induced star formation/positive feedback, give qual- 
itatively similar results. Induced star formation neither leads to 
a fundamentally different response nor causes a global confla- 
gration. Variations in the IMF or introduction of a time-variable 
IMF similarly do not change that picture. In all cases, the vari- 
ations in the SFR can be appreciable, amounting to peak values 
up to 10 times the average value. 

A second type of burst occurs for low gas densities and 
temperatures beyond 10 4 K. This mode is not related to the dy- 
namical evolution, but an instability operating when the cooling 
function drops with increasing temperature. The latter works 
for some temperature regimes (related to H and He line emis- 
sion) above 10 4 K, which can only be reached when there is a 
very low gas density. 

Since the first type of bursts is governed by the evolution of 
the dynamical state of the galaxy, it also seems plausible to dis- 
cuss the results of our simplified models for isolated galaxies 
in other contexts, too. Virial oscillations following structural 
changes induced by galaxy interactions would also lead to an 
enhancement of the SFR, thus providing an additional mecha- 
nism that shapes starbursts. Alternatively, if an interaction does 
not result in a substantial mass redistribution, self-regulation 
might prevent a starburst. 
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